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We introduce a model able to reproduce the main features of 1// noise: hyper-universality (the 
power-law exponents are independent on the dimension of the system; we show here results in 
d — 1,2) and apparent lack of a low- frequency cutoff in the power spectrum. Essential ingredients 
of this model are an activation-deactivation process and dissipation. 
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The voltage drop V on a resistor of resistance R 
through which a current / is flowing obeys the well known 
ohmic law V = RI. Yet, when we look carefully, we 
discover that such a voltage is not perfectly constant 
through time. Indeed there are noise fluctuations around 
V. The spectral density of these fluctuations clearly 
shows a 1// behavior on many decades in the frequency 
domain. This is a well known example of 1// noise, one 
of the most common and widespread features in nature. 
It appears in a variety of systems ranging from the light of 
quasars [jl| to water flows in rivers [pi, music and speech 
[B|, and the already mentioned electrical measurements 
[I,PJ. Despite its ubiquity and universality, a clear and 
simple explanation for such a behavior is still lacking. 
Indeed, it is possible to find in the literature some ad 
hoc formulas and theories, but most of them are based 
on unverified assumptions, or they catch a glimpse of the 
physics only of some particular system, therefore missing 
to address the widespread occurrence of the phenomenon 



In the search for a universal mechanism of 1// noise, 
Bak, Tang and Wiesenfeld (BTW) proposed the new con- 
cept of self-organized critical (SOC) systems J(|: These 
are systems driven by their own dynamics to a state char- 
acterized by power-law time and space correlations, and 
therefore also by power-law (1// Q ) power spectra. Yet, 
a number of features of SOC systems do not show agree- 
ment with the features of 1/f noise: the exponent a is 
seldom close to 1, and it depends strongly on the dimen- 
sionality of the system (at least below the upper criti- 
cal dimension, which is in general high moreover, 
in SOC systems power-law time correlations are always 
found in the presence of power-law (long range) space 
correlations, for which there is no evidence in most sys- 
tems exhibiting 1/f noise [ fio| |. 

In this Letter, we propose a simple model, inspired 
by a SOC model originally introduced by one of us [ 5tJ , 
able to implement some of the current and most accepted 
ideas on 1/f noise, and to show a clear 1/f behavior 
independent on the dimension of the system (therefore 
"hyper-universal" ) . 

The basic model is a continuous version of the BTW 
sandpile. Given a lattice, to every site i is associated a 



continuous variable Xi (representing, say, energy). The 
basic time step of the dynamics consists in changing the 
value of an energy Xi of a positive random quantity e 
(taken from some probability distribution P(e)) 



Xi(t + l)=Xi(t)+e(t) 



(1) 



with the e(t) variables uncorrelated in time. Whenever 
this addition step makes an energy Xi greater than a cer- 
tain value x c , then the quantity Xi is redistributed to the 
2d nearest neighbors of site i 



Xj (t, T + 1) = Xj (t, t) + 



Xjjt, t) 
2d 



(2) 



and the energy Xi is reset to 0. The time variable r is 
used to describe the redistribution process, which is con- 
sidered to be much faster than the process of addition 
(Q). It is possible that this redistribution drives some 
other energies to exceed x c , triggering new redistribu- 
tions. This process (an avalanche in the jargon) goes on 
as long as there are no more energies greater than x c . 
Then a new quantity is added, as in (|l|), and time t in 
increased by 1. It is important to remind the presence 
of two different time-scales: a slow one, corresponding 
to the addition of energy, and a fast one, correspond- 
ing to the redistribution of energies which are above x c . 
The statistical properties of this model after a transient 
time are very interesting. The distribution of the en- 
ergies on the lattice clearly shows a quantization, with 
peaks around 0, x c /2d, 2x c /2d,..., (2d — l)x c /2d. More- 
over the distribution of avalanches with respect to their 
duration (measured on the internal time variable t) obeys 
a power-law. In general it is possible to show that there 
are long-range (power-law) space and time correlations. 
Indeed, this is a SOC model. A further quantity that 
is interesting to look at is the total lattice energy con- 
tent X(t) = Y^,i %i(t)- It represents a signal whose power 
spectrum also obeys some power-law in the frequency do- 
main. 

This model (as well as the original BTW sandpile) is 
mainly an activation/deactivation process, which is be- 
lieved to be one of the main features relevant for the 
description of 1/f noise [fl2|| . Yet, as pointed out above, 
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the exponents depend on the dimensionality of the lat- 
tice, and are never close to 1 [Jul . Therefore the model, as 
it stands, is not a good candidate to describe 1// noise. 

New ingredients need to be added to the model: the 
energy is added only on one side of the lattice, defin- 
ing implicitly a preferred propagation direction for the 
energy; the second ingredient is dissipation. During re- 
distribution, we added some dissipation in the form 



Xj(t, T + 1) = Xj(t, t) + 



2d 



"(I 



(3) 



With this new rule avalanches cannot establish anymore 
long-range correlations throughout the system and are 
not anymore power-law distributed. In a word, dissipa- 
tion destroys the self-organized criticality of the system. 
Yet, some features, such as the quantization of the energy 
levels, survive. 

In our implementation of the model, we inject energy 
on one side of the lattice according to (0), and let it 
propagate through the lattice following (||). We compute 
the power spectrum of X(t\ finding a clear 1/ / behavior 
both in 1 and 2d (see FigJ|) for at least three decades. 
This is a signature of the desired (and observed in nature) 
"hyper-universality" . 

Of course the details of the implementation are rel- 
evant up to some level: We take the added random 
energy e from a uniform distribution in [0, e max \, with 
trnax 3> x c . Indeed, if e rnax < x c then, on the average, 
it takes some steps before the energy goes above thresh- 
old. Since (e(t)e(t')) = 6(t — t'), this implies that short 
time fluctuations of X(t) are uncorrelated and the result- 
ing power spectrum has a flat tail for high frequencies. 
Moreover, we find that for larger systems the frequency 
range where 1// behavior emerges is broader; yet, due 
to the non criticality of the model, the transient time 
to go to stationarity grows fast with the system sizes, 
forbidding us to explore systems larger than 128 x 128 
lattice sites in two dimensions. Also, dissipation can- 
not be too small: indeed, when a becomes very small, 
then we approach the SOC system, which is character- 
ized by different exponents, and crossover effects emerge. 
In Fig.[j] we show the d — 2 power-spectrum with differ- 
ent dissipation regimes. In the absence of dissipation, the 
power-spectrum is essentially flat, with only a small re- 
gion of power-law behavior l/f . A very small dissipa- 
tion (a = 0.0003) clearly gives an intermediate behavior, 
with a high frequency l/f region, and a low frequency 
flat one. Of course, the larger the dissipation, the larger 
the energy injection must be in order to acivate all the 
lattice. 

We believe nonetheless that the way crossover emerges 
is model-dependent. Indeed we believe that the relevant 
features of the model arc nonlincarity (in this case, acti- 
vation/deactivation of sites) and dissipation. 

Additionally, we investigate the power spectrum 
S(f, x) of the energy at site x. We find that, for large x, 



S(f, x) has a scaling form 

S(f,x) = e Sx h(fe 5x ) 



(4) 



From this scaling form we can infer that there is a 
characteristic time T(x) ~ e Sx associated with a site at 
distance x from the origin of the lattice; we can build 
an intuitive picture of this characteristic time thinking 
that in order for the energy to propagate from site x 
to site x + 1, it has to overcome some barrier, with a 
characteristic time to overcome it taken as an Arrhenius 
law e s . Then, in order to propagate from the origin down 
to site x, the characteristic time becomes, roughly, of 
the order of e Sx . As an alternative explanation, due to 
dissipation, energy has a probability to propagate to a 
depth x which is exponentially decreasing with x, hence 
an exponential characteristic time associated with x. 

Due to dissipation, there are no long range correla- 
tions in the system. Therefore, as a first approximation, 
the energies in different sites are uncorrelated. The total 
power spectrum can therefore be written as 

S(/) = X>(/,aO~ [ L e 5x h(fe s *)dx 

fe SL 

= Tfl dyHy) (5) 



(indeed the power spectrum of uncorrelated signals is just 
the sum of the power spectra of the signals, a signature 
of linear superposition). 

We see therefore that the l/f behavior of the power 
spectrum emerges as the superposition of local power 
spectra that have nothing to do with l/f noise. The 



lower cutoff frequency f c 



-SL 



vanishes extremely 



fast in the thermodynamic limit, accounting for the ob- 
served experimental absence of a lower cutoff (whose 
presence is necessary to have a finite power associated 
with the signal). Such a superposition mechanism to 
obtain a 1// power spectrum is strongly reminiscent of 
the McWhorter model Up: S(f) = f S(f, f c )P(f c )df c 
with P(f c ) ~ 1/fc the distribution of the frequencies 
f c and S(f, f c ) ~ f c /(f 2 + fc)- In our model we have 
dfc/fc = dx, accounting for the correct distribution of 
characteristic frequencies, but we have no explicit form 
for S(f, f c ) (although also in our case the large / behav- 
ior is l/f 2 , see Fig||. The relation f c = T~ x = e~ Sx 
can have many different underlying physical origins, as 
diverse as tunnelling between different traps, jumps be- 
tween metastable states distributed in space or dissi- 
pation in an activation/deactivation process (as in the 
present realization) . Such a variety of mechanisms giving 
rise to the good scaling functions (and many other may 
be concieved) strongly points to the observed widespread 
occurrence of l/f noise in nature. 

Actually, dissipation associated with SOC was already 



considered in 14 la] and more recently in uM. There 
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non-trivial power-laws in the power spectrum where 
found, strongly dependent on the dissipation coefficient. 
Moreover dissipation and driving where chosen in such a 
way that the SOC behavior was not destroyed. Actually, 
many experiments show that there is indeed universality, 
that the power spectrum is close to 1/f, and that the 
long range time and space correlations typical of SOC 
systems are absent JTofl . All these ingredients, on the 
contrary, are present in our model, where changing the 
parameter values gives rise to crossover effects but not to 
non-universality. 

Our model shows moreover the relevance of a propa- 
gation direction (feeding energy from one side of the sys- 
tem, and extracting it from the other) , and the 1/f power 
spectrum turns out to be the superposition of power spec- 
tra that are far from 1/f. 

In conclusion, we have introduced dissipation in a 
well-known SOC model: As a consequence, the critical 
behavior of the system is destroyed (but not its self- 
organization properties, such as the quantization of the 
energy levels). The resulting power-spectrum has a clean 
1/f behavior both in 1 and 2 spatial dimensions: We be- 
lieve that our 1 and 2 dimensional results hint toward the 
independence of such behavior on the actual dimension- 
ality of the system, and to the desired hyper- universality. 
We also unveiled that the origin of such a behavior has 
to be found in the superposition of power-spectra with 
characteristic frequencies f c suitably distributed in space. 
Such a distribution is not an input in the model, but 
emerges due to the directedness properties of our model. 
The desired distribution of characteristic times typical 
of the McWhorter model emerges spontaneously in our 
model. One can concieve many other situations where 
the distribution of the characteristic frequencies is in- 
stead given, and the system only has to cope with it. 
As a consequence, the present model not only provides a 
simple and hyper-universal explanation of 1/f noise, but 
is also suggestive of a wide variety of microscopic physical 
mechanisms able to give 1/f noise. 
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FIG. 1. Log-log plot of the power spectrum of X(t) in d — 1 
and d = 2. The straight lines are 1/f power-laws drawn for 
reference. The system size used for the d = 1 simulations is 
L — 100, with €max = 10 and a = 0.03; d — 2 simulations 
are performed with L = 100, e max = 1 (energy is added on 
all the sites xo,i, for a maximum possible energy injection of 
100) and a = 0.01. For comparison, we also add the d — 2 
power-spectrum with a much smaller dissipation (a — 0.0003, 
but with the same size and energy injection regime as before) 
that shows crossover between clean 1/f behavior, and the be- 
havior in the absence of dissipation (lowest curve) that clearly 
shows no sign of 1/f behavior. 
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FIG. 2. Log-log plot of the power spectra S(f, x) for 
x = 25, 35, 45 on a d — 2 lattice of L — 100, e max = 10 
and a = 0.01. The large / behavior has a characteristic 
slope l// 2 . The collapse in the inset is obtained plotting 
e~ Sx S(f, x) vs. fe Sx with 8 = 1.9. Indeed the universality of 
the scaling function (4) emerges. 
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